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CN I We present a study of the nearby Seyfert galaxy NGC 1068 using mid- and far- 

infrared data acquired with the IRAC, IRS, and MIPS instruments aboard the 
Spitzer Space Telescope. The images show extensive 8 /xm and 24 /xm emission 
coinciding with star formation in the inner spiral approximately 15"(1 kpc) from 
the nucleus, and a bright complex of star formation ~47" (3 kpc) SW of the 
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^ I nucleus. The brightest 8 /im PAH emission regions coincide remarkably well 

with knots observed in an Ha image. Strong PAH features at 6.2, 7.7, 8.6, and 
1L3 yum are detected in IRS spectra measured at numerous locations inside, 
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within, and outside the inner spiral. The IRAC colors and IRS spectra of these 
regions rule out dust heated by the AGN as the primary emission source; the 
SEDs are dominated by starlight and PAH emission. The equivalent widths 
and flux ratios of the PAH features in the inner spiral are generally consistent 
with conditions in a typical spiral galaxy ISM. Interior to the inner spiral, the 
influence of the AGN on the ISM is evident via PAH flux ratios indicative of a 
higher ionization parameter and a significantly smaller mean equivalent width 
than observed in the inner spiral. The brightest 8 and 24 /im emission peaks 
in the disk of the galaxy, even at distances beyond the inner spiral, are located 
within the ionization cones traced by [O III]/H/3, and they are also remarkably 
well aligned with the axis of the radio jets. Although it is possible that radiation 
from the AGN may directly enhance PAH excitation or trigger the formation of 
OB stars that subsequently excite PAH emission at these locations in the inner 
spiral, the orientation of collimated radiation from the AGN and star formation 
knots in the inner spiral could be coincidental. The brightest PAH and 24 /xm 
emitting regions are also located precisely where two spiral arms of molecular gas 
emerge from the ends of the inner stellar bar; this is consistent with kinematic 
models that predict maxima in the accumulation and compression of the ISM 
where gas gets trapped within the inner Lindblad resonance of a large stellar bar 
that contains a smaller, weaker bar. 

Subject headings: galaxies: Seyfert — galaxies: individual (NGC 1068) — in- 
frared: galaxies 



1. Introduction 



Polycyclic aromatic hydrocarbons (PAHs) appear to be ubiquitous in the ISM of massive 
disk galaxies. PAHs emit copiously near regions of active star formation where molecules are 
heated by ultraviolet radiation and are usually associated with the photo-dissociation regions 
in molecular clouds. The infrared spect ra of active galaxie s generally lack the PAH emission 
l ines c haracteristic of starburst galaxies (iRoche et al.lll99ll ). iFreudling. Siebenmorgen. fc Haas 
(120031 ) reported the detection of hot dust and weak PA Hs in type 1 AGN, and cooler dust and 



strong PAH s in type 2 AGN, simi lar to the results of IClavel et al.l (120001 ). iBuchanan et al. 



(120061 ) ■ and lAppleton et al.l (120071 ). These results are interpreted as an indication that the 
high energy photons from the AGN destroy PAH molecules exposed to extreme UV radiation 
from the nucleus, with a ny observed PAH emission coining from star forming regions shielde d 
by intervening material (jVoitlll992bl : lAllain et al.lll996l : lMaloneylll999l : IFreudling et al.ll2003l ) . 
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The Spitzer Space Telescope provides the opportunity to investigate the distribution of PAH 
emission in a variety of galaxies through its correlation with the zone of influence of the 
AGN and/or star- forming regions. 



NGC 1068 is a prototypical nearby Sey fert 2 galaxy (lAntonucci fc Milled Il985l ). and 
was the first AGN to be detected (jFathlll909t ). With an infrared luminosity Lir = 10^^'^ Lq, 
NGC 1068 is a luminous infrared galaxy (LIRG) and is a member of the IRAS Revised 
Bright Galaxy Sample (RBGS; Sanders et al. 2003) consisting of extragalactic objects with 
60/im flux greater than 5.24 Jy. In addition to the AGN, the galaxy contains a bright 
inner spiral with active star formation approximat ely 1 kpc from the nucleus; this structure 



has often been referred to a s a "starburst rin g " (|Mvers &: Scovilld 119871 : ISchinnerer et al. 
2000; Spinoglio et al. 2005 ). Telesco &: Decher ( 19881 ) estimated that star for mation would 



Yuan&Kuo 



consu me the available gas within 0.5 Gyr at current star formation rates 
(119981 ) showed that the inner spiral resides at the outer Lindblad re sonance. The inner 
spiral arms begin at the ends of a stellar bar JSchinnerer et al. 2000 ). at the location of 
peak collisional energy dissipation (jCombes fc GerinI Il985l ). Fabry- Perot imaging of the 
[OIII]A5007 line shows clearly defined ionization cones from the center of NGC 1068, inclined 
such that the disk of t he g alaxy is illuminated o u t to r ~ 11 kpc ( Fig. [H adapted from 
Bruhweiler eTaP J200lh and lVeilleux et all J2003h; lUnger et al.lll992l ; IVeilleux et aLll2003h . 



The [Nil] A6548, 6583 hne profile data of ICecil et al.l (Il990l ) showed that the kinematics of 
the ionized gas are well fit by a model in which the ionization cones intersect the plane 
of the galaxy, rather than illuminating halo gas seen in projection onto the disk. The 
brightest region of [OIII] emission corre sponds with the region of peak radio co ntinuum 
flux, extending ~ 10" NE of the nucleus JCondon et al.lll99ol : [Morganti et al.lll993h . A knot 
of [OIII] emission is seen at the distance of the inner spiral, with further [OIII] emission 
extending several kpc further out into the disk of the galaxy. Kinematics show that the gas 
in the ionization cones is ionized in situ by the central engine, and is not ejecta from the 
nuclear regi on. Spectral line rat ios also confirm that the AGN is the ionizing source even at 
r ~ 11 kpc (IVeilleux et al.ll2003l ). though the inner spiral itself shows line ratios typical of star 
form ing regions. The o p ening angle of the iori ization cones is 40°, centered at position angle 
31° (lUnger et al.lll992l ). IVeilleux et al.l (120031 ) showed that the [OIII] emission is biconical, 
though mu c h fain ter to the SW than the NE. This is consistent with the conclusion of 
Cecil et al.l (Il990l ) that the ionization cones are inclined with respect to the plane of the 
disk, illuminating the near side of the disk to the NE and the far side of the disk to the 
SW. iBruhweiler et al.l (12001! ) showed that the [0III]/H/3 ratio both interior and exterior to 
the inner spiral along the cones is co nsistent with photo ionization from the AGN without 
significant contribution from shocks. lUsero et al.l (120041 ) detected the inner spiral in SiO, 
which is typically interpreted as a signature of shocks. Infrared emission peaks have been 
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Bland-Hawthorn et al.l (119971 ) argued that the IR emission results from dust grains near the 
surface of mole c ular c louds heated by the X-ray and UV flux from the AGN. Using ISO, 
Le Floc'h et al.l (120011 ) detected the peaks in the 7.7 /xm PAH band, indicating that at least 
some of the emission is from PAHs rather than large dust grains. 

Is it possible that the inner spiral is not illuminated by the AGN, but simply appears 
to be illuminated due to projection effects? If the inner spiral lies in the plane of the 
disk, the AGN ionization cones would have to be inclined well out of the plane of the disk 
and have a small opening angle so as not to intersect it. As noted previously, the [Nil] 
and [QUI] data sets provide strong evidence that the ionization cones intersect the disk. 
Packham et al.l (119971 ) used near-IR polarimetry to produce a three dimensional model of 
the electron scattering cone, concluding that it is inclined by 43° to the plane of the sky, 
corresponding to an inclination of < 10° to the plane of the galaxy. With an opening half- 
angle of 40°, it clearly intersects the disk. The remaining possibility is that the inner spiral 
is inclined relat i ve to the disk such that the inner spiral is not illuminated by the AGN. The 
Packham et al.l (119971 ) model requires an inclination of at least 30° for the inner spiral to lie 
outside of the ionization cone. Such an inclination is both physically implausible and would 
have been easily detected in kinematic data (e. g. CO or Ha). We therefore conclude that 
the AGN illuminates the inner spiral in the region of the PAH knots. The fact that optical 
line ratios characteristic of AGN ionization are seen both immediately interior and exterior 
to the inner spiral indicates that material in the inner spiral cannot be completely shielded 
from the AGN by intervening material. 

In this paper we present the results from analysis of new infrared obs ervations of 



NGC 1068 from the Spitzer S pace Telesco pe. Observations with the IRAC (IFazio et al. 



2004f ). IRS (iHouck et al.ll2004j ). and MIPS (iRieke et al.ll2004j ) instruments are discussed in 
§ 2. Analysis of the data is presented in § 3, and conclusions are given in § 4. 



2. Observations and Data Analysis 
2.1. Imaging 

NGC 1068 was observed with IRAC by the GTO program of Fazio et al. (PID #32). 
MIPS observations were carried out as part of the Great Observatories All-sky LIRG Survey 
(GOALS, GO-1 RBGS LIRG/ULIRG program of Mazzarella et al; PID #3672). The general 
data reduction procedures used on the IRAC and MIPS data were described in detail in 
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Mazzarella et al.l (120071 ) . Here we discuss the special processing required for NGC 1068. For 
reference, the bright star forming regions are at a radius of 15" (1 kpc). 

Muxbleed saturation artifacts (see the IRAC Data Handboolj^ for details) in the 5.8 /im 
and 8 /im images presented challenges in the nuclear region of NGC 1068 within a radius 
of 15" (1 kpc). A rough correction for this problem consisted of measuring each artifact in 
a suitable aperture and subtracting off the flux density of an identical aperture measured 
at the same galaxy surface brightness level in a region outside the artifact. This procedure 
resulted in a correction (reduction) of the flux density inside a 15" aperture by 2.74 Jy (16%) 
at 8.0 /im and 0.41 Jy (5%) at 5.8 fim. 

Since muxbleed is a result of saturation, an additional correction was needed to account 
for the lost flux density. Such saturation corrections were estimated by comparison with 
available data in the literat ure. The flux density measured within a 6" dia meter aperture 



was 3. 2 ± 0.2 Jy at 5.0 /xm (IRieke fc Lowlll972al ) and 25 ± 3 Jy at 10.0 /xm (IRieke fc Low 



1972 bl ). A simple linear interpolation from these values to the IRAC wavelengths predicts 



6.2 ± 0.5 Jy and 16 ± 2 Jy at 5.8 and 8.0 /xm, respe ctively. The blue en d of the transmissio n 



window of the fllter used by iRieke fc Lowl (Il972bl ) extends to 7.6 /im (ILow fc Riekelll97ll ) 



so the 10 /im measurement includes the majority of the flux density from the 7.7 /im PAH 
feature. The Spitzer flux densities measured in a 6" diameter aperture centered on the 
nucleus of NGC 1068 are 7.9 Jy at 5.8 /im and 7.0 Jy at 8.0 /im. (We note that the flrst 
muxbleed artifact in the 8 /xm PSF is located at a radius of about 6", which is well outside 
this comparison region with a 6" diameter.) The 5.8 /im value measured from our corrected 
IRAC data agrees with the estimate based on the Rieke & Low measurements within the 
uncertainties of our correction method. However, since our 8 /im IRAC estimate is 9 Jy 
below the estimate based on interpolation of the Rieke & Low data, we have added 9.0 Jy 
to the 8 /xm IRAC aperture measurements for NGC 1068 to account for the flux density lost 
to saturation. 

Photometry of the inner spiral was measured by taking the difference between measure- 
ments of 12" and 18" radius (800-1200 pc). Extended source corrections were applied to 
each aperture measurement according to the equation derived by T. Jarrettj^. The emission 
peaks on the inner spiral were measured using small elliptical apertures. The wings of the 
PSF of the bright AGN would bias the observed flux densities if the emission peaks were 
measured relative to the sky background, so instead the peaks were measured relative to 
a local background elsewhere on the inner spiral. In the case of the NE knot, the local 



^http: / / ssc.spitzer.caltech.edu / irac /dh/ 

^http: / /spider. ipac.caltech.edu/staff/jarrett/irac/calibration/ext_apercorr. html 
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background measurements were chosen to also lie along diffraction spikes at 3.6 and 4.5/im. 

The MIPS 24 /im image of NGC 1068 is dominated by a bright unresolved central point 
source (the AGN) which saturated within radius ~ 13". To measure the spatially resolved 
features we subtracted a PSF model of the central point source. The IDL package IDPsi was 
used to centroid the galaxy and PSF images; then the scale factor factor of the PSF image was 
interactively adjusted to best subtract off the point source component of the galaxy image. 
As recommended by the MIPS Data HandboolJ^, an object of similar color was chosen to 
create the PSF image. The ULIRG IRAS F05189-2524 served as a satisfactory PSF for this 
purpose, as it is both very bright and unresolved by Spitzer at 24/im. IRAS F05189-2524 is 
also a Seyfert 2 galaxy, and is the best match to the IRAS 60 /im/25 /xm color of NGC 1068, 
with a color of 3.82 compared to 2.24 for the latter. The only galaxies in the GOALS sample 
with colors more closely matching NGC 1068 were spatially resolved in the MIPS 24 /im 
images. Contours of the PSF-subtracted 24 /im image of NGC 1068 are overlaid on the 
8 /im image in Fig. [2l The elongation of the masked region on the 8 /im image is due to the 
muxbleed residuals discussed above. 

The best fit PSF for NGC 1068 was also used to estimate the total 24 /im flux density, 
which could not be measured directly due to saturation. The point-source component of the 
total flux density was estimated by taking the flux density of IRAS F05189-2524 (3.04 Jy) 
and multiplying by the best fitting PSF scale factor (25), resulting in an estimate of 75.90 Jy. 
In addition, the extended 24 /im emission from NGC 1068 was measured directly from the 
PSF-subtracted image. Combining the 0.5 Jy extended component with the point-source 
component yields an estimated total 24 /xm flux density of NGC 1068 (before aperture and 
color corrections) of 76.4 Jy. The uncertainty in the empirically determined PSF scale factor 
(25.0) is about 20%, and this dominates the uncertainty in this 24 /zm total flux density 
measurement. 

Table 1 presents th e Spitzer flux density m easurements and estimated uncertainties. 



As described in detail in iMazzarella et al.l (120071 ). MIPS color corrections have been applied 



based on the IRAS and MIPS SEDs. With aperture and color corrections, the total estimated 
24 /im flux density using the Spitzer data is 79.6 Jy. Comparison with the IRAS 25 /im 



measurements of 87.6 ±0.2 Jy fISanders et al.ll2003l ) indicates that this 24 /im estimate from 



Spitzer is low by about 9%. This is consistent with the systematic mean offset of 20% 
(with a Icr scatter of 15%) observed between the 24 /im Spitz er measurements and 2 5 /im 



IRAS measurements across our entire sample of 203 systems (IMazzarella et al.l 120071 ). We 



■^http: / /mips. as. arizona.edu/mipspage/IDP3/idp3example. html 
^http:/ /ssc. spitzer. caltech.edu/mips/dli/ 
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are therefore confident in the PSF fitting procedure, from which it can be stated that more 
than 99% (75.9/76.4) of the 24 /xm emission from NGC 1068 originates within a radius of 
^ 30" (the first Airy ring). Approximately 80% of the 24 /im emission originates within the 
saturated region r < 13". NGC 1068 also saturated the 160 fim array and a flux density 
estimate could not be recovered. 



2.2. Spectroscopy 

The nucleus of the galaxy was observed at high spectral resolution by the GTO program 
of Houck et al. (PID #14), and mapped at low resolution by the GO-1 program of Gallimore 
et al. (PID #3269). The Short-Low (SL) IRS maps consisted of 13 half slit-width (1.8 arcsec) 
steps for the first and second order slits. The maps were centered on the nucleus. There 
were no dedicated off (sky) observations, but a local background was subtracted by using 
the non-primary slit, which is centered approximately 80 arcsec from the nucleus, for both 
the SLl and SL2 maps. The IRAC 8 fim surface brightness at the locations selected for 
the background are typically 10~^ - 10~^ of the level at the nucleus. CUBISMlf] was used to 
perform the local background subtraction and form the SL spectral cubes. From these cubes 
we were able to identify regions over which to extract ID spectra, and we used CUBISM to 
perform this extraction in matched apertures in SLl and SL2. In addition, we used CUBISM 
to generate a pure PAH map over the area covered by the IRS slits. Aperture locations and 
the PAH map are shown in Fi g. [3l Fluxes and e quivalent widths of the low resolution data 



were measured using PAHFIT (jSmith et al.ll2007l ) . The regions closest to the nucleus suffered 



from saturation and did not yield usable spectra; see Fig. [3l lower panel. Table 2 presents 
the equivalent widths and line fluxes of the emission features measured in each aperture. 

The Short-High and Long- High IRS observations do not suffer from saturation, but only 
cover a single position at the center of the galaxy. No background subtraction was performed 
due to the lack of a background slit observation. Spectra wer e extracted using SP IClj^, and 



lines were measured on the extracted spectra using SMART flHigdon et al.ll2004l ). Our high 
resolution spectral line measurements are presented in Table 3. 



^http://turtle. as. arizona.edu/jdsmith/cubism.php 
^http://ssc.spitzer.caltech.cdu/postbcd/spice.html 
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2.3. Separating Non-Stellar from Stellar Emission in Seyfert Galaxies 



Although IRS spectra provide the ideal measure of the presence and strength of PAH 
emission, the available data cover only a portion of the inner spiral. IRAC and MIPS imaging 
is required to investigate the region outside of the spectral map. There are several possible 
sources for the IR emission from Seyfert galaxies. The power-law continuum from an AGN 
can contribute to all IRAC wavelengths. The 3.6 /im and 4.5 /im data are dominated by 
the Rayleigh- Jeans tail of stellar photospheric emission. The 5.8 /im channel covers the 
wavelength range which includes the 6.2 fim PAH line, while the 8 fim channel includes the 
7.7 fim and 8.6 fim PAH emission features. Hot dust can also be a significant contribution 
at all IRAC wavelengths. Stellar continuum emission is ins ignificant compared to PAH 
emission in the IRAC 8 /im band in nearly all nearby LIRGs (IMazzarella et al.l 120071 ). The 
stellar continuum flux density at 8 u m is estimated to be 19-23% of the 3.6 /im flux density 
(jPahre et al.ll2004j : iDale et al.ll2005l ): based on this we estimate that stars contribute only a 
few percent (for example, NGC 1068 1 kpc aperture: 2.8 x 0.23/17 = 0.038) of the 8 um flux 
densities of this galaxy. To test for spatial variations between the starlight and PAH emission 
in NGC 1068, a scaled 3.6 /im image was subtracted from the 8 /im image. The ratio of 
this nonstellar emission image to the 8 /im image is constant with little variation across the 
galaxy (excluding the central region which saturated at 8 /im, as noted previously), showing 
that starlight accounts for ~ 4 % of the 8 fim emission throughout the disk of NGC 1068. 

Determining the presence of PAHs using IRAC and MIPS photometry requires the 
measurement and modeling of the contribution of AGN continuum and hot dust at 8 /im. 
Several color-color diagnostics have been used to d istinguish AGN-d ominated systems and 



starburst-dominated sy stems in the literature (e. g. iStern et al.ll2005l ). The GOALS sample 
( IMazzarella et al.l 120071 ) is shown in the [3.6]-[4.5] vs. [4.5]-[8] color-color space in Fig. |H 
Also shown are the colors of hot dust and AGN power laws. Since the vertical axis is the 
ratio of two channels normally dominated by starlight, galaxies without significant AGN or 
hot dust lie near zero. The horizontal axis separates galaxies with little or no PAH emission 
(left) from those dominated by PAH emission (right). The [4.5]-[8] color measures the jump 
from the Rayleigh- Jeans tail of starlight at 4.5 /im (plus AGN continuum, if present) to the 
bright PAH emission at 8 /im. The better-known Stern et al. color-color space uses [5.8]-[8] 
on the horizontal axis instead of [4.5]- [8]. Since both the 5.8 and 8 /im channels include PAH 
bands, [4.5]-[8] is preferred due to the clearer physical interpretation. 

Combinati ons of IRAC i mages can be used to create spatial maps of different components 
of IR emission. iPahre et al.l (120041 ) used 3.6 and 4.5 /im data to estimate the contribution of 
starlight, which was then scaled down in proportion to a Rayleigh- Jeans tail and subtracted 
from the 8 /im image to create a map of nonstellar (typically PAH) emission. As noted 
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above, this technique has httle effect on LIRG and ULIRG data such as those examined here 
since the 8 /im images are dominated by nonstellar emission. A similar method can be used 
to identify nonstellar emission at 4.5 fim. No PAH emission lines are present at such short 
wavelengths, so the nonstellar emission sources are direct radiation from the AGN and hot 
dust only. Making an assumption about the form of the nonstellar emission, this technique 
can be used to quantitatively estimate the emission from stars and the assumed nonstellar 
source at 3.6 an d 4.5 /zm. The rat io of stellar emission at 3.6 /im to that at 4.5 /im, f^, 
is known; using iPahre et al.l (120041 ) = 1.7. A similar ratio, /pi is easily calculated for 
any given power law slope a. For typical values 0.84 < a < 1.4, 0.83 < /pi < 0.73. 
Image subtraction is then performed, with the power law component of the 4.5 /im image 
given by /pi(4.5) = ^ where /(A) is the IRAC image at a given wavelength. 



Similar equations are easily derived for the stellar components of each image, or the power 
law component of the 3.6 /xm image. Note that if the nonstellar emission arises from hot 
dust this method will only reveal the locations of such emission and one can simply use 
7(4.5) — 7(3 .6) //h. without attempt ing to account for the nonstellar contribution to the 
3.6 /xm flux (lEngelbracht et al.ll2005l ). A quantitative flux density estimate would require an 
assumed dust temperature from which /dust could be calculated, replacing /pi in the above 
prescription. In either case, the assumed nonstellar emission can also be extended to 8 /im 
to better constrain the PAH emission. 



3. Discussion 



3.1. Spitzer Data 



The IRAC colors of the region within 1 kpc of the nucleus and the total emission from 
the galaxy place N GC 1068 in the AGN region of the IRAC color-color diagram defined by 
Stern et al.l (120051). Emission at IRAC wavelengths is highly concentrated, with 70-80% of 
the flux density coming from the central kpc (Table 1). 

High resolution spectra of the nucleus show strong AGN signatures in a number of di- 
agnostic line ratios. Line fluxes and equivalent widths are presented in Table 3. 
The [O IV]25.9 /im/[Ne II]12.8/im line ratio is u sed as a diagnostic for the fractional c ontri- 
bution of the AGN to the luminosity of a galaxy (IGenzel et al.lll998l : lArmus et al.ll2007l ). The 
NGC 1068 nucleus has a ratio of 4.2 as measured from the IRS spectra, compared to 2.7 as 
measured from the ISO data of lSturm et al.l (120021 ). [Ne V]14.3 /im/[Ne II]12.8/im is a similar 
diagnostic ratio; the nucleus of NGC 1068 has a ratio of 2.1, compared to the IS O ratio of 1.4. 



Altho ugh these AGN diagnostic ratios are higher than the ISO measurements of lSturm et al. 



( 120021 ). this is not surprising since the IRS slits are smaller than the corresponding ISO aper- 
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ture by roughly a factor of two. With less light from the disk of the galaxy, the line ratios 
might be expected to show higher ionization. Using the [Ne V]14.3 /im/[Ne III] 15.6 fim and 
[Ne II]12.8/xm/[Ne III] 15.6 una line ratios measured using the high resolution spectrum, the 
diagnostic diagram of IVoitI (Il992al ) indicates an ionization parameter log U ~ —2 and a 
power law ' input spectrum of index a ~ —1.3. Compared to nearby AGNs, starbursts, and 
ULIRGs (lArmus et al.l 120071 ). the line ratios in the nucleus of NGC 1068 are in excellent 
agreement with the AGN sample, consistent with a 100% AGN contribution to the observed 
mid-IR nuclear spectrum. 

The IRAC images show bright 8 fim peaks (also detected in the other IRAC channels) 
along the N E and SW arcs at PA ^ 45° in the inner spiral (Fig. [2]) at the locations first 
detected by iTelesco fc Decher Jl988h . Of the total 8 /im flux density within the annulus 
800 < R < 1200 pc from the nucleus, 52% is emitted by these two knots. The surface 
brightness of the knots (400-500 MJy/sr) is more than twice that seen elsewhere in the 
inner spiral (~ 200 MJy/sr). Compared to the central kpc of the galaxy, these knots are 
bluer by more than half a magnitude in [3.6]-[4.5] and redder by 0.55-0.85 mag in [4.5]-[8]. 
Since the wavelength range of the 8 nm channel includes both the 7.7 /xm and 8.6 /xm PAH 
features, the latter color strongly suggests that the knots are dominated by PAH emission. 

The low resolution IRS data confirm the predominance of PAH emission in the knots 
observed in the IRAC 8/xm image. The 6.2 /im, 7.7 fim, and 11.2 fim PAH features, plus the 
PAH12.7/im+[Ne II]12.8/im blend, are clearly visible throughout the inner spiral. Figure O 
shows representative spectra, including the inner spiral, the region interior to the inner spiral, 
and the region outside of the inner spiral. Fluxes and equivalent widths of each measured 
emission feature are presented in Table 2. The PAH features are strongest in the inner 
spiral (apertures 1-5, 9-11; top panel of Fig. [5]). However, such features are also present in 
the innermost usable apertures (apertures 16, 18, & 19; middle panel of Fig. E]). Although 
the PAH equivalent widths in these innermost apertures are small due to the bright hot 
dust continuum near the AGN, the fluxes of the PAH lines are several times as bright as in 
aperture 12 (outside the inner spiral) and typically half as bright as in the inner spiral. These 
innermost apertures also show the strongest [S IV] 10.51 /im emission, while [Ar II] 6.98 fim 
and H2 9.66 fim emission are strongest on the inner spiral. 



For comparison, the average starburst galaxy spectrum from iBrandl et al.l (120061 ) was 
also analyzed using PAHFIT. Line fluxes and equivalent widths measured using PAHFIT 
cannot be directly compared with measurements made using other techniques b ecause PAH 



FIT t akes into account the flux in the extended wings of the spectral lines; see ISmith et al. 



J2007h for details. PAH equivalent widths are 2.36 /xm, 7.15 fim, 1.09 /im, 1.36 /xm, and 
0.97 fim (6.2, 7.7, 8.6, 11.3, and 12.7 /xm features, respectively). Though the inner spiral is 
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the source of the strongest PAH emission in NGC 1068, this emission is significantly weaker 
than that seen in starburst galaxies. 

With the detection of bright PAH emission lines throughout the inner spiral, the model 



of iBland-Hawthorn fc VoitI (119931 ) in which the IR peaks on the inner spiral originate from 
dust grains heated by the AGN can be ruled out. A composite model in which the IR peaks 
are due to a combination of PAH emission plus continuum from hot dust heated by the 
AGN can also be strongly constrained. As shown in Fig. [21 apertures 1 and 9 both lie within 
the ionization cones and are the closest apertures to the PAH knots on the inner spiral. If 
hot dust made a significant contribution to the flux in these apertures, the PAH equivalent 
widths would be lower than in apertures on regions of the inner spiral falling outside the 
ionization cones (apertures 3-5, 10, and 11). Instead, the PAH equivalent widths measured 
in apertures 1 and 9 are as high or higher than those measured elsewhere on the inner spiral. 
Any contribution of continuum flux from hot dust heated by the AGN must be too small 
to significantly affect the PAH equivalent widths. We estimate the amount of additional 
continuum flux within the ionization cones compared to the rest of the inner spiral by 
comparing the flux ratios and equivalent width ratios of the PAH features. PAH fluxes in 
apertures 1 and 9 are brighter than the average values in the rest of the inner spiral by factors 
of 20-80%; equivalent widths are brighter by factors of 10-100%. The signature of a greater 
contribution from the continuum is a larger increase in flux than in equivalent width. This 
is seen in aperture 9, with a flux increase of 31% and equivalent width increase of 15% in 
the 11.3/im and 12.7/im features that would be most affected by hot dust continuum. These 
ratios indicate a continuum increase of 1.31/1.15 = 14% compared to the average value in 
the rest of the inner spiral. However, in aperture 1 the increase in PAH equivalent width is 
greater than or equal to the increase in PAH flux for every feature, indicating a continuum 
level equal to or less than that seen in the rest of the inner spiral. Comparing the aperture 1 
and 9 spectra to the average spectrum of the inner spiral apertures outside the ionization 
cones confirms the conclusion that little or no excess continuum is seen within the ionization 
cones. Figure [6] shows the three spectra normalized to a common flux density at 13 fim. The 
apertures within the ionization cones clearly show brighter PAH features than the average 
spectrum, while the continuum shows very little variation between the three spectra. Warm 
dust visible at 24 fim (Fig. [2]) is present at the locations of the peaks on the inner spiral, 
though it is unclear whether the AGN or local star formation activity is the primary heating 



source. 



Ratios be tween the 6.2 / xm, 7 .7 /zm, and 11.3 PAH features are shown in Fig. [71 Lines 



adapted from iDraine fc Lil (120011 ) show the locations of PAH emission from a cold neutral 
medium (CNM) model, a warm ionized medium (WIM) model, and a photodissociation 
region (PDR) model. Also shown is a portion of the boundary of the region of PAH line 



- 12 - 



ratios which can be reproduced by the composite ISM model of iDraine fc Lil (120011 ) . With 
the exception of aperture 10, the PAH measurements of the inner spiral all lie above the 
boundary, within the region consistent with a composite model. Aperture 5 is the outlier 
in the direction of a colder environment and smaller PAH molecules. The H2 9.66/im line 
is nearly a factor of two brighter in this aperture than in the other apertures. Interestingly, 
aperture 18 — only 480 pc from the nucleus — also lies in this cold ISM region of the 
diagram. It is also interesting that aperture 10 on the inner spiral and aperture 12 outside 
the inner spiral both lie in the region requiring more ionizing radiation than the composite 
ISM model, with similar PAH ratios to the majority of the apertures interior to the inner 
spiral. Since apertures 10 and 12 are both well outside the [OIII] ionization cone, ionization 
from the AG N is unlikely. Ape r ture 1 coincides with the location of the peak CO emission 
in the map of ISchinnerer et al. J2000h . suggesting the formation of massive stars likely 
ionization source. The general picture presented by the ratios of PAH features is that PAH 
molecules in the inner spiral are smaller and lie in a cooler, less ionized environment than 
PAHs closer to the nucleus. 

The CNM-like PAH ratios of aperture 18 suggest that radiation from the AGN may 
be absorbed within a few hundred parsecs along the NE ionization cone. The radio jet of 
Condon et al.l (Il990l ) ends abruptly near the same location. However, optical line ratios also 
show that radiation from the AGN is ionizing gas up to 11 kpc from the nucleus within the 
NE cone jVeilleux et J1 1200,4 § 3.2). 



Although the 24 fim MIPS data are saturated (r < 13", 900 pc), the innermost available 
isophotes are elongated along the same NE-SW axis corresponding to the arcs seen at other 
wavelengths (Fig. [2]). In addition to the PAH knots, a weaker feature is visible in both the 
IRAC and MIPS 24 fim images, corresponding to the NE tip of the southern arm of the inner 
spiral (22" from the nucleus at PA = 64°). Two additional peaks are seen ~ 45" (3 kpc) to 
the SW of the nucleus at PA = 214°. The nonstellar emission map (Fig. [S]) shows that hot 
dust is present in these locations. The subtraction of the nuclear PSF (right panel) suggests 
that there is little hot dust emission in the inner ~ 20" outside of the nucleus; widespread 
heating of dust to ~ 1000 K can be ruled out. It is curious that every peak of IR emission 
in t he disk of NGC 1068, in cluding the three such locations on the inner spiral as defined by 
the ISchinnerer et al.l (120001 ) CO map, lie within the ionization cones of the AGN. 



3.2. Comparison to Other Wavelengths 



The radio continuum emission ( Condon et al. 1119901) comes to an abrupt end at a distance 
of 500-700 pc from the nucleus, roughly halfway to the inner spiral, as shown in Fig. O The 
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Sum knots on the inner spiral lie ~ 5" (300 pc) farther along the radio axis. iCecil et al. 
(Il990l ) found that the high velocity nuclear outflow ends abruptly at the NE radio lobe and 
the SW radio hotspot (4.5" from the nucleus). 



The CO observations of lSchinnerer et al.l (120001 ) (Fig. [9]) show a coherent spiral structure. 
The 8 /zm emission peaks coincide with regions of bright CO emission, though the converse 
is not true. The two brightest PAH knots are seen at precisely the locations where the radio 
jet axis intersects the inner spiral (Fig. [H]). The 200-500 pc distance between the ends of the 
radio jets and the inner spiral indicate that the AGN is unlikely to trigger shock-induced 
star formation, however. 



The optical emission line imaging of IVeilleux et al.l (120031 ) provides important insight 
into the origin of the observed PAH knots. The PAH knots coincide with strong Ha emission 
(Fig. [ini left panel) and with small [O lll]/H/? ratios (0.2-1.7; Fig. ^ right panel). The 
8 /xm and Ha images show strikingly similar morphologies, tracing active star forming regions 
across a wide range in wavelength. The [O lll]/H/3 ratio is an indicator of mean ionization 
and temperature; high ratio s (> 3, as a rule of thumb) ar e typi cal of AGN, while low er ratios 
are typical of H ll regions (IVeilleux fc Osterbrockl 119871 ) . As IVeilleux et al.l (120031 ) pointed 
out, high [O lll]/H/3 ratios are seen out to i? ~ 11 kpc in NGC 1068, most visibly in the NE 
ionization cone. Figure [TO] clearly shows that both the starburst knots and the three fainter 
IR emission peaks have optical line ratios indi cative of star fo r matio n. Although the Ha 
and [O lll]/H/5 data confirm the conclusion of iLe Floc'h et al.l (120011 ) that the PAH knots 
are due to star formation in those locations on the inner spiral, it is still surprising that 
bright PAH emission should be present from regions that are being illuminated by ionizing 
radiation from the AGN. 

While the galaxy disk is illuminated by the AGN, it is unclear how opaque the PAH 
emitti ng regions are to radia tion from the AGN (e.g. see Fig. [TTl left panel). The GALEX 
data (iGil de Paz et al.l l2006l ) provide little guidance on this issue. Figure [11] (right panel) 
is a GALEX FUV image of NGC 1068 shown with contours of both the Spitzer 8/im and 
PSF-subtracted 24/im emission overlaid. As is apparent in the figure, the regions showing 
strong PAH and warm dust emission are clearly UV luminous, indicating that (i) the PAH 
emitting regions are optically thin to UV photons and thus unshielded by the AGN, that (ii) 
the UV emanates from the optically thin outer envelope of a much more embedded, optically 
thick starburst responsible for producing the observed PAH emission, or (iii) that the UV 
emission is AGN light scattered by optically thick clouds. The latter possibility is unlikely 
— a simple assumption of a 100 pc spherical cloud with an albedo of one would reflect a 
flux ~ 10^ fainter than the flux of the AGN. Comparing measurements within 6" (GALEX 
FWHM) radius apertures centered on the nucleus and on each PAH knot, the GALEX data 
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show that the PAH knots are ~ 10% as bright as the nucleus in the FUV and therefore 
scattering from the AGN is unhkely to be a significant contribution. Thus, the UV emission 
is associated with stars, but cannot be used to rule out whether the PAH regions are shielded 
from the AGN. 

Further interpretation is complicated by the near-coincidence between the axis of the 
AGN jet and the pos ition angle of the inner stellar bar. Intense star formation is common 
at the ends of a bar (jCombes fc GerinI 119851 ). so the AGN may have no effect on the PAH 
knots at all. It is also possible that the AGN is either triggering increased star formation or 
directly exciting PAH emission. Either the high energy photons capable of destroying PAHs 
are absorbed near the AGN while lower energy photons are unobstructed over the entire 
~ 11 kpc of the illuminated disk, or the high energy photon flux within 1 kpc produces 
a PAH destruction rate less than or equal to the PAH formation rate in the star-forming 
regions. For details on the p arameters governing the formati on and destruction o f PAH 
molecules, see iMaloneyl ( 19991 ). The radiative transfer model of ISiebenmorgen et al.l ( 20041 ) 
supports the latter hypothesis, finding that even in the optically th in regime PA Hs survive 
exposure to a 10^^ Lq AGN at distances of 100 pc or more. However, IVoitI (1l992bl ) calculated 
that PAHs could only survive direct exposure to a 2.5 x 10^° Lq AGN at distances > 14 kpc. 



4. Conclusions 

We have investigated the influence of the AGN on PAH emission in NGC 1068, a LIRG 
with a Seyfert 2 nucleus. Bright PAH knots are detected in the inner spiral in NGC 1068. 
These knots lie directly along the axis of the ionization cone, which extends well beyond the 
inner spiral in the disk of the galax y. The knots also coinci d e wit h the ends of the inner 



stellar bar. Contrary to the model of iBland-Hawthorn fc VoitI (Il993l ). the spectra and colors 



of the PAH knots are consistent with starlight and PAH emission, with very little if any hot 
dust. Three other star forming regions are also detected: the NE tip of the southern arm of 
the inner spiral, and two regions 2.5-3 kpc to the SW. All three regions lie within the AGN 
ionization cones and show PAH emission at 8 /im and a slight excess of 4.5 fim emission 
relative to starlight, suggesting dust heated by young stars. IRS spectra near the PAH knots 
indicate that PAHs are the dominant energy source, and show that PAH emission is present 
within a few hundred parsecs of the AGN. PAH line ratios on the inner spiral are generally 
consistent with a cool, neutral ISM, while PAH line ratios interior to the inner spiral generally 
indicate a more ionized medium. These data show that the AGN is not destroying all PAH 
molecules within the ionization cone. The enhanced PAH emission within the ionization 
cones indicates that the AGN is either triggering increased star formation, directly exciting 
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the PAHs, or has no effect on the dust properties in the ionization cones. 
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Table 1. Spitzer Flux Density Measurements 



IRAC 



MIPS 



Feature 


Aperture Radius 


3.6 /im 


4.5 fim 


5.8 fiia 


8 fim 


24 /im 70 fMia 






axcRcc (kpc) 


•ly 


■h' 


■h' 


■h' 


•ly .ly 


NGC 1068 


Total 


90" (6 kpc) 


3.8 


5.1 


13 


23 


80 180 


NGC 1068 


1 kpc 


15" 


2.8 


4.3 


9.0 


17 




NGC 1068 


NE PAH knot 


4" (0.26 kpc)* 


0.011 


0.006 


0.05 


0.11 




NGC 1068 


SW PAH knot 


5" (0.33 kpc)=i 


0.024 


0.020 


0.11 


0.26 




NGC 1068 


Inner Spiral 


12"-18" (0.8-1.2 kpc) 


0.29 


0.20 


0.80 


1.85 





^Effective radius of custom aperture. Flux densities for the PAH knots are relative to the local background elsewhere 
in the inner spiral. 



Note. — The approximate uncertainties are indicated by the number of significant figures; specifically, 10% at 
IRAC wavelengths, and 20% at 70 /um. See text for a discussion of the 24 uncertainties. 
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Table 2. Emission Features (Low Resolution) 



Aperture 


r (kpc) 


PAH 6.2 


[Ar II] 6.98 


PAH 7.7 


PAH 8.6 


H2 9.66 


[S IV] 10.51 


PAH 11.3 


PAH 12.7 


[Ne II] 12.8 


































Inner Spiral 












1 


0.95 


328 i 6 


12 it 2 


1020 ± 20 


175 ± 3 


9.0 ± 0.8 


5.1 ± 0.6 


285 ± 5 


170 ± 2 


35.7 ± 0.4 


(4.3", 13.5") 




2.04 


0.042 


3.91 


0.65 


0.066 


0.032 


1.08 


0.60 


0.19 




0.88 


250 ± 4 


9 ± 2 


780 ± 10 


118 ± 3 


7±1 


6.8 ± 0.6 


228 ± 2 


124 ± 2 


22.0 ± 0.4 


(0.4", 13.2") 




1.71 


0.066 


3.14 


0.43 


0.043 


0.039 


0.80 


0.42 


0.11 




0.87 


268 ± 8 


8 lb 3 


820 ± 30 


130 ± 10 


11 ± 3 


5 ± 2 


226 ± 8 


94 ± 8 


22.9 ± 0.9 


( — 5.4", 11.9") 




0.55 


0.024 


1.73 


0.32 


0.026 


0.029 


0.75 


0.29 


0.11 




0.90 


236 ± 4 


9 ± 2 


770 ± 10 


127 ± 4 


10 ± 1 


3.9 ± 0.6 


245 ± 3 


112 ± 2 


24.2 ± 0.5 


( — 10.1", 9.0") 




1.22 


0.058 


2.85 


0.49 


0.088 


0.028 


0.96 


0.36 


0.12 




0.98 


221 ± 4 


6 ± 3 


560 ± 20 


153 ± 4 


19 ± 2 


5 ± 1 


275 ± 6 


100 ± 2 


20.8 ± 0.4 


( — 14.0", 4.5") 




1.03 


0.031 


1.69 


0.57 


0.21 


0.042 


1.20 


0.31 


0.10 




0.98 


268 ± 3 


13 ± 2 


950 ± 10 


143 ± 3 


4.5 ± 0.8 


7.7 ± 0.6 


287 ± 2 


134 ± 2 


28.5 ± 0.4 


(-2.5", -14.4") 




2.04 


0.089 


3.50 


0.48 


0.030 


0.043 


0.92 


0.39 


0.12 


10 


1.02 


125 ± 3 




644 ± 9 


82 ± 2 




4.2 ± 0.3 


143 ± 1 


55 ± 1 


10.7 ± 0.3 


(4.3", -14.7") 




0.93 




4.84 


0.47 




0.021 


0.42 


0.13 


0.038 


11 


1.02 


148 ± 3 


5 ± 1 


570 ± 10 


79 ± 2 


2.9 ± 0.6 


5.6 ± 0.4 


199 ± 1 


127 ± 1 


13.4 ± 0.3 


(14.8", -3.8") 




0.67 


0.032 


2.18 


0.30 


0.021 


0.035 


0.78 


0.48 


0.075 










Inside the Inner Spiral 












0.68 


133 ± 3 


6 ± 2 


630 ± 20 


105 ± 3 


2.0 ± 0.5 


4.2 ± 0.4 


155 ± 2 


86 ± 2 


15.0 ± 0.4 


(-10.1", -1.8") 




0.24 


0.018 


1.29 


0.25 


0.008 


0.019 


0.50 


0.30 


0.079 




0.65 


179 ± 6 


4 ± 2 


830 ± 20 


133 ± 4 


8.1 ± 0.8 


12.0 ± 0.5 


244 ± 3 


100 ± 3 


23.7 ± 0.4 


(-7.2", -6.5") 




0.22 


0.008 


1.12 


0.22 


0.030 


0.039 


0.46 


0.15 


0.053 




0.74 


126 ± 4 


1.7 ± 0.9 


660 ± 20 


137 ± 3 


4.1 ± 0.5 


7.2 ± 0.4 


223 ± 2 


117 ± 2 


19.1 ± 0.3 


(-5.4", -9.7") 




0.20 


0.004 


1.20 


0.29 


0.016 


0.029 


0.60 


0.29 


0.069 


13 


0.68 


157 ± 4 


5 ± 1 


670 ± 20 


120 ± 3 


4.3 ± 0.4 


7.2 ± 0.4 


171 ± 2 


101 ± 2 


17.1 ± 0.3 


(2.2", 10. 0") 




0.37 


0.017 


1.56 


0.30 


0.018 


0.031 


0.48 


0.26 


0.064 


14 


0.66 


146 ± 5 


5 ± 2 


640 ± 20 


103 ± 3 


2.8 ± 0.4 


6.3 ± 0.4 


150 ± 2 


91 ± 3 


15.5 ± 0.4 


(-2.2", 9.7") 




0.32 


0.019 


1.56 


0.26 


0.011 


0.024 


0.39 


0.25 


0.063 


15 


0.64 


141 ± 5 


5 ± 2 


620 ± 20 


104 ± 4 


3.3 ± 0.5 


5.4 ± 0.4 


125 ± 2 


71 ± 3 


16 ± 6 


(-6.8", 6.8") 




0.30 


0.016 


1.47 


0.26 


0.012 


0.021 


0.33 


0.19 


0.066 


17 


0.78 


95 ± 3 


1.8 ± 0.7 


480 ± 10 


83 ± 2 


1.6 ± 0.4 


4.2 ± 0.3 


116 ± 1 


54 ± 2 


10.6 ± 0.3 


(11.5", -1.6") 




0.23 


0.006 


1.19 


0.22 


0.007 


0.019 


0.39 


0.19 


0.056 










Apertures Nearest the Nucleus 










16 


0.54 


99 ± 5 


2 ± 1 


630 ± 20 


41 ± 4 


7.4 ± 0.3 


28.6 ± 0.7 


187 ± 2 


48 ± 3 


18.7 ± 0.7 


(-2.2", -7.8") 




0.11 


0.003 


0.64 


0.04 


0.010 


0.031 


0.11 


0.02 


0.012 


18 


0.48 


165 ± 4 


12 ± 1 


690 ± 10 


52 ± 3 


2.5 ± 0.6 


32.5 ± 0.4 


296 ± 2 


306 ± 3 


34.0 ± 0.8 


(6.1", 3.7") 




0.24 


0.024 


0.84 


0.05 


0.003 


0.030 


0.16 


0.14 


0.023 


19 


0.54 


119 ± 4 


6 ± 2 


600 ± 20 


88 ±4 


2.7 ±0.5 


5.3 ± 0.5 


133 ± 2 


92 ± 2 


11.8 ± 0.6 


(-7.9", 1.7") 




0.18 


0.015 


0.98 


0.18 


0.011 


0.022 


0.34 


0.19 


0.035 










Outside the Inner Spiral 










12 


1.37 


35 ± 2 




173 ± 7 


26 ± 2 


0.7 ±0.3 




45.8 ± 0.8 


16 ± 1 


2.3 ± 0.2 


(17.3", -11.1") 




0.90 




2.22 


0.26 


0.010 




0.42 


0.15 


0.031 



Note. — Fluxes are in units of 10 "^""^ W cm ^, with the equivalent width in microns below. 
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Table 3. Nuclear Emission Features (High Resolution) 



Line 


Flux 


Equivalent width 




10-21 W cm-2 


/Lim 


[S IV] 10.51 


460 ± 60 


0.006 


[Nc II] 12.81 


520 ± 60 


0.008 


[Nc V] 14.3 


1110 ± 60 


0.017 


[Nc III] 15.55 


2100 ± 200 


0.024 


[S III] 18.71 


220 ± 20 


0.004 


[Ne V] 24.31 


800 ± 100 


0.023 


[O IV] 25.89 


2200 ± 200 


0.077 


[S 111] 33.48 


420 ± 60 


0.025 


[Si II] 34.81 


620 ± 20 


0.034 


[Ne III] 36.0 


130 ± 20 


0.006 



1. — Diagram of NGC 1068 ionization cone. Left panel: view in the plane of the disk, 



from Bruhweiler et al 



(120011) ■ Righ t panel: RO continuum image with ionization cones illustrated 



(adapted from IVeilleux et alJ (120031 ) ). Large tick marks are at 1' intervals. 
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Fig. 2. — The 8 fim image of NGC 1068 (grayscale), with contours of the PSF-subtracted 24 nm 
image. The brightest contour levels are shown in light gray for clarity. The saturated data in the 
nucleus have been excluded from both images. Residual banding is visible on the 8 //m image as 
the streak running SE-NW. The ring of spots is also an artifact of the 8 /xm PSF. The 24 /xm data 
show many of the features seen at 8 /xm, including the two knots on the inner spiral, the NE tip of 
the southern arm of the inner spiral 1.5 kpc (23") ENE of the nucleus, and the peak 3 kpc (45") 
SW of the nucleus. 
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Fig. 3. — Extraction apertures for the Sliort-Low IRS spectra are overlaid on the 5.8 IRAC 
mosaic (top, green numbered boxes). Also shown are the area covered by the Short-Low map 
(blue box) and the position of the Short-High slit (yellow box). Banding and PSF artifacts are as 
described in Fig. [2j North is up and east is to the left. The map of the 6.2 fim PAH feature is 
shown in the bottom panel. The center of the galaxy has been masked out. 
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Fig. 4.— IRAC colors of the RBGS LIRGs and ULIRGs from Mazzarella et alJ ((20071). NGC 1068 
is marked (1 kpc radius and total apertures; red and black points respectively), as are the colors 
of the two bright knots on the inner spiral and the colors of the galaxy disk (r > 1 kpc) and the 
inner spiral. Small red points indicate the central 1 kpc radius aperture measurements of the rest 
of the GOALS sample, while black points are total magnitudes. Unit vectors indicating the change 
in color from the central aperture to the total aperture are shown; note that there is not a one 
to one correspondence between central and total aperture measurements due to a combination of 
resolution and multiple interacting galaxies, respectively. The median change in color from the 
central 1 kpc aperture to the total aperture is shown by the arrow in the upper right. Colors 
are expressed in the Vega magnitude system, so by definition a pure starlight SED would appear 
near zero in both color axes. Also shown are the colors of a pu r e AG N power-la w spectrum (open 
square s, from lower left to upper right: a = 0.84, IClavel et al.l (120001 ): a = 1.4, iNeugebauer et al 
jiasa)), and hot dust (blue points for 1000 K and 900 K). 
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Fig. 5. — IRS Short-Low spectra of NGC 1068. The top panel shows the average of eight apertures 
on the inner spiral, the middle panel shows the average of the three apertures nearest the nucleus, 
and the bottom panel shows aperture 12, outside the inner spiral. 
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Fig. 6. — Low resolution spectra of the inner spiral within the ionization cones (aperture 1, blue; 
aperture 9, red) are compared with the average spectrum of the inner spiral outside the ionization 
cones (black). The three spectra have been normalized to a common flux density at 13 /xm. No 
difference is found between the continuum flux density within the ionization cones compared to 
outside the ionization cones. 
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Fig. 7. — Relative strengths of three PAH emission features. Apertures along the inner spiral are 
shown as open circles. Apertures inside the inner spiral are shown as solid circles The aperture 
exterior to the inner spiral is shown as an open triangle. Outliers have been labeled by aperture 
number. Un certainties ar e ~ .01 in each PAH flux ratio, and are not shown. The lines are 
adapted from lDraine &: Lil (|200ll ). The solid line corresponds to the CNM model, the dashed line 
corresponds to the WIM model, and the dotted line corresponds to the PDR model. For all three 
models, PAH molecule size increases to the left. The red dotted line shows part of the boundary 
of the region above which a simple ISM model can reproduce observed PAH features. Data below 
this line require greater ionization. 
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Fig. 8. — Left: Contours of an estimate of the nonstellar emission overlaid on 4.5 //m grayscale 
image of NGC 1068. As expected, the AGN dominates the nonstellar emission, but the PAH 
emission regions outside the inner spiral (SW and ENE of the nucleus; marked) are also visible. 
The method used to estimate the nonstellar component of the 4.5 flux density is described 
in § 2; Q = 1 was used as the power law slope. Right: The 4.5/im nonstellar emission from the 
point-like nucleus has been removed by PSF subtraction. The circle marks the location of the inner 
spiral 1 kpc (15") from the nucleus. This image shows that there is very little detectable emission 
from hot dust outside of the nucleus. 
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Fig. 9. — The 8 //m image of the center of NGC 1068, with 1.49 GHz radio continu um contours from 
Condon et alj (1990) (red) and ^^CO (1-0) contours from lSchinnerer et alj ( 2000l ) (white) overlaid. 



Note that the peaks of 8 /xm emission he at the locations where the radio jet axis intersects the 
inner spiral. 
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Fig. 10. — Left: Ha image of NGC 1068 f rom IVeilleux et al.l (120031 ) overlaid with 8 ^m contours 



Right: [O iii]/H/? image of NGC 1068 from lVeiheux et al.l (|20o3 7 overlaid with 8 fim contours. The 



orientation of the ionization cones is shown, and the PAH knots on the inner spiral are labeled. 
The [O ili]/H/3 ratio is low throughout the inner spiral (the light gray regions on the inner spiral 
have ratios ~ 1). Elsewhere in the ionization cone, particularly to the NE, the [O iii]/H/3 ratio 
is several times as large (the dark gray regions have ratios ~ 7), indicating that the AGN is the 
ionizing source. 
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Fig. 11. — Left: Diagram showing the optically thick (top) and optically thin (bottom) scenarios 
described in the text. The diagram is not to scale. Blue lines represent UV photons from the AGN 
or massiv e stars; red lines represe nt PAH emission. Right: The GALEX FUV (1516A; grayscale) 
data from I Gil de Paz et al.l ()2006l ) is compared with Spitzer 8/xm (blue) and PSF-subtracted 24;um 
(red) data. The regions of the inner spiral showing strong PAH and warm dust emission are also 
bright in the UV, indicating that these star forming regions are not optically thick to UV photons. 
The slight offset (~ 3" NNE) between the GALEX and MIPS emission is real; the vahdity of 
the GALEX WCS solution was confirmed by plotting the positions of USNO catalog stars on the 
GALEX image. 



